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Abstract 

We demonstrate how one should transform correctly quasi-isotropic coordinates to Weyl- 
Papapetrou coordinates in order to compare the metric around a rotating star that has been 
constructed numerically in the former coordinates with an axially symmetric stationary metric 
that is given through an analytical form in the latter coordinates. Since a stationary metric as- 
sociated with an isolated object that is built numerically partly refers to a non- vacuum solution 
(interior of the star) the transformation of its coordinates to Weyl-Papapetrou coordinates, which 
are usually used to describe vacuum axisymmetric and stationary solutions of Einstein equations, 
is not straightforward in the non- vacuum region. If this point is not taken into consideration, one 
may end up to erroneous conclusions about how well a specific analytical metric matches the metric 
around the star, due to fallacious coordinate transformations. 

PACS numbers: 95.30. Sf, 04.25.D-, 04.40.Dg 



1 



I. INTRODUCTION 



Recently there has been quite a few attempts to describe the geometry around an astro- 
physical object, such as a rotating neutron star, a strange star, or a black hole surrounded 
by an accretion disk, through various types of analytical solutions of the vacuum Einstein 
equations [5|. Numerous people have produced a lot of families of such exact 

solutions during the last decades following various generating techniques. The metrics that 
correspond to these families of solutions are usually parameterized by a few parameters. 
These parameters could then be used to fit the specific characteristics of each type of central 
object by appropriate tuning of the parameter values. Thus if the matching between a large 
range of metrics that are constructed numerically and a family of parametrized analytical 
metrics is quite acceptable, the space-time neighborhood of a multi-parametric central com- 
pact object could be represented quite accurately by this family of exact solutions. One 
could then use these metrics to explore analytically the behavior of orbits around such a 
central object (compute the frequencies related with these orbits, find out the innermost 
circular radius, etc.), or the other way around: probe the physical characteristics of the 
central body itself by exploiting the characteristics of the gravitational waves produced by 
small objects orbiting around the central one. 

Fortunately, nowadays, there is a large variety of analytical solutions of vacuum Einstein 
equations, which could be used as candidate metrics to describe well the exterior space-time 

n 

of axisymmetric astrophysical objects. Ernst [6] formulated the Einstein equations in the 
case of axisymmetric stationary space-times long time ago, while Manko et al. and Sibgat- 
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121 ] have used various analytical methods to produce such space-times 



parameterized by various parameters that have a different physical context depending on the 
type of each solution. Also, Neugebauer [3] have constructed a specific axisymmetric solu- 
tion analogous to the well known Schwarzschild and Kerr black hole solutions: it describes a 
a rotating thin disk of dust. Among all these available solutions, one has to choose a specific 
type of solution that relates better to the particular astrophysical object, depending on the 



specific physical characteristics one expects 
On the other hand various groups (see 



rom such an object. 



sec 



14J ] , and for an extended list of numerical schemes 



151]). that have expertise in building relativistic models of astrophysical objects with 



adjustable physical characteristics, can construct the metric inside and outside such objects 
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by solving numerically the full Einstein equations in stationary cases. Their numerical codes 
generate metrics in tabulated form with numerical values that correspond to the metric 
components at the grid points that have been assumed in their numerical scheme. 

In order to compare an analytical solution with a metric that has been constructed nu- 
merically, one should make sure that the transformation of the coordinates of the metrics to 
each other, if not the same, is absolutely faithful. More specifically a problem arises when 
one attempts to transform the quasi-isotropic coordinates, that are usually used to describe 
a metric that has been constructed numerically, both inside and outside the astrophysical 
object to Weyl-Papapetrou coordinates, that are usually used in the analytical expressions 
of the available stationary axisymmetric metrics which will be compared with the numerical 
one. Although the transformation of coordinates is straightforward in the vacuum region, 
the same type of transformation leads to erroneous coordinates when it is used in the mat- 
ter region. This inconsistence is expected since an axisymmetric metric in quasi-isotropic 
coordinates is described by four metric functions, while in the Weyl-Papapetrou coordinates 
there are only three independent functions present in the usual stationary axisymmetric met- 
ric (actually not all of them are independent, since in both cases one of the metric functions 
is uniquely obtained, apart from a constant, from the rest functions). The difference comes 
actually from the fact that the former one is used for an axisymmetric stationary solution 
of the general Einstein equations, while the latter one is used for axisymmetric stationary 
vacuum solutions of Einstein equations. Of course when one deals with vacuum solutions 
in quasi-isotropic coordinates the four metric functions are interrelated to each other by an 
extra constraint. On the other hand if one insists on using the Weyl-Papapetrou coordi- 
nates p and z inside the matter, one could not anymore use them on equal basis as in the 
usual vacuum axisymmetric stationary metric; instead one more function A(p, z) should be 
introduced to describe the induced metric of the two-dimensional surface spanned be p and 
z, namely 



In this paper we suggest that in order to translate the former to the latter coordi- 
nates, the corresponding integration path that is usually used to compute the z— coordinate 
(c.f. Sec. IHIj) should avoid entering the matter region. Moreover, since the exterior region 
is the one we are interested to in order to exploit its characteristics to probe the source of 




(1) 
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the gravitational field, we are not really interested to know the actual coordinates inside the 
matter. 

The rest of the paper is organized as follows: In Sec. [IT] we show the relation between the 
two sets of coordinates and the recipe to compute the metric components in Weyl-Papapetrou 
coordinates from the metric components in quasi-isotropic coordinates. In Sec. II III we argue 
that the best path of integration to compute the z— coordinate is to follow a meridian 
(r = const), starting from an equatorial point just outside the star's surface, up to whatever 
angle 9 and then move along the radial coordinate r either inwards up to the surface of the 
star or outwards to infinity. We end up this section by a practical formula for computing 
the z— coordinate. Furthermore we use this formula to obtain the exact relation between the 
z— coordinate and the isotropic r— coordinate in the Schwarzschild metric as a demonstration 
of of the proposed scheme. In Sec. IIVI we give an estimate of the errors that arise from 
numerical integration of z, and the corresponding errors that are induced in the metric 
components. Once again we use the example of Schwarzschild space-time to measure these 
numerical errors since in this very case we know exactly the z-coordinate and we can compare 
it with the value of z obtained numerically. We argue that the errors computed for the 
Schwarzschild case are of the same order of magnitude (~ 1CT 6 — 1CT 7 in g tt ) as for any 
neutron star model obtained by the numerical code of Stergioulas, even for the most rapidly 
rotating ones, if the simple trapezoid rule is used in numerical integration. Thus we conclude 
that any relative difference between a numerical and an analytical metric of order higher 
than 10~ 6 should be attributed to a real non- matching of the metrics. In the last section 
we summarize our conclusions and show how much improved is the comparison between the 
numerical and analytical metrics studied by Berti and Stergioulas if the transformation 
of the coordinates is done according to our proposed scheme. 

All physical quantities used in this paper are in geometrized units (G = c = 1). 

II. COORDINATE TRANSFORMATION IN VACUUM. 

The line element of an axisymmetric and stationary space-time in quasi-isotropic coordi- 
nates assumes the following form: 

ds 2 = -e 2u dt 2 + e 2 ^ (d<f) - uodtf + e 2 » (dr 2 + r 2 d6 2 ) , (2) 
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where v, ip, u, p are functions of r and 9 alone. An alternative way to write this metric is 
by replacing the ip function by a new function B through 

= Be~ u rsin6 (3) 

upon which the line element transforms to: 

ds 2 = -e 2u dt 2 + B 2 e~ 2v r 2 sin 2 9 (d<f> - udt) 2 + e 2 <+ (dr 2 + r 2 d9 2 ) . (4) 

We note once again that the metric written in any of the above forms has full freedom to 
describe any stationary axisymmetric solution of the full Einstein equations. 

On the other hand the Weyl-Papapetrou coordinates are very good to describe any sta- 
tionary axisymmetric solution of the vacuum Einstein equations. The line element in these 
coordinates is 

ds 2 = -f (dt - wd(j)f + f- 1 [e 27 (dp 2 + dz 2 ) + p 2 d<j) 2 ] , (5) 

where now the three functions /, w, 7 are functions of p, z alone. 

In order to transform the quasi-isotropic coordinates to Weyl-Papapetrou coordinates 



(c.f. 



one first defines the cylindrical coordinates 

zu = r sin 6, ( = rcos9. (6) 

Then the Einstein field equations in vacuum (R\ + Rt = 0) imply (c.f. [if]]) 

d 2 (^B) d 2 {wB) 

dw 2 dc, 2 K) 

Thus one could use a new coordinate 

p = wB = e (l/+ ^ (8) 

instead, which satisfies the two-dimensional Laplace equation in the (r — 9) surface. One 
could then define a harmonic function that is conjugate to p, that is 

z = z(w,0, (9) 

which satisfies the Cauchy-Riemann conditions 



At this point we should note that the above construction of conjugate coordinates p, z was 
feasible only in the vacuum region. In the interior of the axisymmetric star, p could still be 
defined as above but it does not anymore satisfy the two-dimensional Laplace equation and 
thus the other coordinate z fails to be constructed as an harmonic conjugate of p. 

Now in the vacuum region of space-time we can integrate the Cauchy-Riemann conditions 
(flQj) with initial value for z 

z(w, C = 0) = 0. (12) 

This corresponds to the equatorial plane of the star, or to be more specific, to that part of the 
equatorial plane that lies outside the star. The corresponding integration yields the value of 
z— coordinate at any point outside the star. Although z = at the equatorial plane inside 
the star as well, we cannot integrate these relations along a path that lies inside the star 
since they do not hold in matter. Thus the path we choose to integrate the Cauchy-Riemann 
relations should lie entirely in the vacuum region up to the final point. 

Besides transforming the (r, 9) — or (w, Q — coordinates to (p, z) coordinates, one has to 
compute the new metric functions (f,w,j) from the old metric functions (y, u),p, B) that 
are supposed to be known at the grid points used in the numerical code that generates them. 
By direct use of the Cauchy-Riemann relations (llOflip . and subsequent substitution of the 
coordinates defined in relation ([6]) we get 

2 



dp 2 + dz 2 



dp 
dw 

dp 
dw 



dp 
dC 

dp 



{dw 2 + d( 2 ) 
(dr 2 + r 2 d6 2 ). 



(13) 



Finally by comparing the two metrics (141151) . and keeping the coordinates (t, 0) the same in 
both metrics, we obtain 
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dp 

dw 



dp 



(14) 
(15) 

(16) 



These new metric functions, that are computed from the old ones, along with the new co- 
ordinates — p that is directly computed from the old coordinates, and z that is computed 
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by integration along paths that lie entirely along vacuum regions — complete the metric 
transformation in Weyl-Papapetrou coordinates. Of course the integration of the Cauchy- 
Riemann relations in order to compute the z— coordinate cannot, in general, be performed 
analytically; one should rely on some numerical scheme to integrate the corresponding rela- 
tions. This technical issue will be addressed in the following Section. 



III. TRANSFORMATION TO WEYL-PAPAPETROU COORDINATES. 

In this section we will present a practical recipe that one could follow to integrate the 
Cauchy-Riemann relations in order to obtain the Weyl-Papapetrou z— coordinate. As men- 
tioned in the previous section, the other conjugate Weyl-Papapetrou coordinate, p, is di- 
rectly obtained from the metric function B and the w— coordinate (cf. Eq. (jSJ)). Through 
the Cauchy-Riemann relations (Eqs. fllOlllip ) it is easy to verify that 

dz dB 

— = cos 9B + sin#— , (17) 
or do 

dz 2 . „ (dB B\ 

Moreover since p = cos# (not to be confused with the corresponding metric function in 
quasi-isotropic coordinates) and x = r/(r + r e ), where r e is the equatorial radius of the star, 
are usually used in numerical schemes to produce numerical values of the metric functions 
in a specific grid that is more uniformly distributed in azimuthal angles and covers better 
the whole space up to infinity, the above derivatives could straightforwardly be transformed 
to 

?1 = r e (x 2 — + B—) . (20) 
dp e \ dx 1 — x ) 

Equipped with these expressions, we may now choose a suitable path to integrate them 
in order to assume the numerical values of the new z— coordinate. In realistic rotating stars 
the equatorial radius of the star r e is the maximum value of r s {6), that is the function that 
describes the shape of the surface of the star. Thus we could simply start from an equatorial 
point {6 = 7r/2) just outside the surface of the star where z = 0, and follow the grid points 
along the meridian x = xq (r = ro — r e ), until we reach the axis of symmetry (9 = 0). 
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Upon reaching whatever intermediate angle 6, we could then move radially (along r, or x) 
— either outwards, or inwards up to the surface — to obtain the numerical value of z at every 
grid point in the vacuum region where the metric is known. We should stress at this point 
that along the meridian of constant x = Xq the optimum integration is achieved by using 
expression ()2Up . while along the radial direction at constant 9 (or /i) the optimum integration 
is achieved by using expression (ITT)) instead of (fl9l) ; this difference in effectiveness of the two 
relations with respect to radial changes of z— coordinate arises because the corresponding 
integrant is well behaved in the former case as r — > oo, in contrast to what's happening in 
the latter case as x — > 1. Thus the z value at the grid point (x, fi) is the outcome of the 
integral 




The above expression simplifies much if the meridian path, followed initially, corresponds to 
exactly the equatorial radius, since then xq — 1/2 (ro = r e ). Practically though, since the 
numerical computation of the derivative dB/dx needs at least two neighboring grid points 
lying in the vacuum region, the meridian path should correspond to the r value of the next 
after the first grid point lying outside the equator of the star. Finally the last integral term 
in the expression above could be omitted if we seek to compute the z values at the grid 
points along the axis of symmetry since there \i — 1. 

Of course one could choose any other path starting from the equator where z = to 
reach the final grid point, but since the integration will be carried numerically it is better 
to choose a path that minimizes the numerical errors. In the next section we will show why 
the path suggested above is expected to be efficient with respect to numerical errors and 
we will give an estimate of the error magnitude. Heuristically, the basic argument in favor 
of this path is the fact that if we follow to move along another meridian x = const which 
is far outside the surface of the star, the error in the numerical computation of dB/dx will 
be much greater, since r(x) has a rapidly increasing derivative as x — > 1. This numerical 
error will then follow as a systematic error in all z values when the integration along x is 
computed next. 

We will end this section by demonstrating this coordinate transformation by a very simple 
example where the integrations could be performed analytically; namely the Schwarzschild 
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metric. In the Appendix one could find the form of the line element of this metric in 
isotropic coordinates, and indirectly in Weyl-Papapetrou coordinates. The B function for 
the Schwarzschild metric is 



Thus, only the dB/dx and B parts will survive in the integrands of Eq. (l2Tj) . For example 
along the positive part of the z-axis (fi = 1) the value of z is easily computed to 



By a few simple substitutions back and forth between the coordinate r and the compactified 
coordinate x we get at the end 



This is the exact relation between the two sets of coordinates (see Appendix). In the next 
section we will use once again the Schwarzschild example to get a first estimate of the errors 
arising from numerical implementation of the recipe described above. 

IV. ESTIMATE OF NUMERICAL ERRORS 

In order to check whether a metric given in analytical form describes faithfully the metric 
of a physically realistic configuration which has been constructed through numerical schemes 
and thus its components are given in tabulated form at the grid points where the metric 
has been computed, one should check how well the two metrics coincide at these specific 
grid points. In the usual case, where the two types of metrics refer to different kind of 
coordinates, one should first transform one set of coordinates to the other set so as to compare 
the two metrics at the same points and then decide about the faithfulness of the specific 
analytical metric. However the transformation of coordinates involve numerical errors since, 
in practice, the integration associated with the computation of the new coordinates will in 
general be performed through numerical integration. Especially if we have the case presented 
in the previous sections where the isotropic coordinates of the numerical metric have to be 
transformed to Weyl-Papapetrou coordinates of the analytical metric, the integration will 
be based on the numerical metric that is given in discrete form and thus the implemented 




(22) 






(24) 
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errors could not be optimized further than a minimum value related to the number of the 
grid points. Therefore one should first have an estimate of these errors before evaluating the 
pure differences between the two metrics. Hence if the differences computed between the 
two metrics are of the order of the numerical errors induced in the metric by the numerical 
transformation of coordinates, we could not assess any countable difference between the two 
metrics. 

An alternative way to tackle the problem of comparing the two types of metric, that 
although it implies numerical errors these could in principle be minimized at will, is to 
transform the coordinates in which the analytical metric is expressed to the other set of 
coordinates. Since the former metric is analytically known, the corresponding integrations 
used to compute the new coordinates could be performed at whatever level of accuracy one 
desires. In that case the difference between the metrics at the same grid points at which 
the numerical metric is known is true and does not correlate at all with the computation 
of the new coordinates. In the case considered in this paper this kind of transformation of 
coordinates is not that simple as the inverse one. Actually there is no direct way to compute 
the isotropic coordinates (r, 6) from the Weyl-Papapetrou coordinates (p, z) as well as the 
metric expressed in the latter ones. Therefore we will resort in the first method, that was 
described in detail in the previous section, to transform coordinates. We will show though 
that the errors induced from numerical integrations, at least for the number of grid points 
used in a specific physical example, is not much larger than the accuracy at which the 
numerical metric is known. This fact suggests that there is no harm in using this method 
to transform coordinates and then compare metrics. 

Let's say then that we want to compare two metrics, a numerical one #^9 ( r >^)> that 
corresponds to a rotating neutron star with specific internal physical characteristics, and an 
analytical one g^J (p, z) which we believe, or simply want to check if, it describes quite well 
the former numerical one. As is shown in the corresponding variables of the two functions 
the two metrics are assumed functions of different coordinates. We decide to compare the 
two metrics along the axis of symmetry of their axially symmetric space-time. Anyway 
the comparison along the axis of symmetry seems to be more demanding for the numerical 
transformation of coordinates since the integration path from the equator to the z-axis 
is longer than to any other angle, and therefore it imposes larger error contribution from 
the first two integrals of Eq. (12ip . On the other hand the last two integral terms in (12ip 
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are competitive to each other for various values of the angle parameter /i. Although we 
cannot draw a general rule about the angle parameter \i at which these two terms assume 
the highest total value, we have seen in practice that the z-axis is really the most heavily 
infected direction from numerical errors in computing the z-coordinate. 

In Figure [TJ we have plotted the logarithmic relative error in numerical computation of 
the z-coordinate in Schwarzschild metric along the z-axis (the metric used corresponds to a 
spherically symmetric star that has the same mass M and equatorial coordinate radius r e as 
the most rapidly rotating neutron star model for which we want to estimate the errors caused 
in g tt ). The number of grid points assumed are the ones used in all numerical metrics that 
lave been constructed for various models of rotating neutron stars by Berti and Stergioulas 
l| and a simple trapezoid method of integration is used. Even though this example could 
not be considered suitable to check the errors in a realistic rotating neutron star case (due to 
spherical symmetry of Schwarzschild metric B is not a function of //, hence the computation 
of the first two integrals in Eq. (pill do not contribute any error) , it gives at least a minimum 
estimate of the order of magnitude of errors. In the spherically symmetric case the errors 
arise from the numerically estimated value of dB/dx in the first integral of Eq. ( I2TI) . and 
the numerical computation of the third integral (the last integral is zero in every case along 
the z-axis). Although the ^-coordinate that we compute numerically gets shifted more and 
more dramatically as z increases due to cumulative errors along the integration path, the 
error induced in the metric itself (e.g. the gu component) does not increases so much with z. 
This is expected since the metric becomes less sensitive to z as we recede from the neutron 
star. Consequently since the metrics themselves are the ones that we want to compare, the 
cumulative error in z at large values of z is not disturbing. In Figure [1] we have plotted 
the corresponding logarithmic relative error in g u , that is due only to erroneous numerical 
integration of z coordinate, along the z-axis. As shown in the plot the relative error does 
not even exceed the ~ 10 -6 , which is actually just about one order of magnitude higher than 
the level of accuracy of the numerical metrics produced by the numerical code of Stergioulas 
[yj]. Thus we conclude that if the size of the rest of errors which are coming from the 
[i— dependent terms of a realistic numerical model do not exceed the errors arising in the 
simple Schwarzschild case there is no need to worry about any numerical errors induced to 
metrics caused by transformation of coordinates. 

In the remaining part of this section we will argue that this is exactly the case with 
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FIG. 1: The plot shows in a log-log plot the relative difference between the exact z value (z^) and 
the z-coordinate produced through numerical integration (z^ N >) at fixed grid points (upper thick 
dashed curve). The model used to estimate the numerical error is a Schwarzschild metric that has 
the same gravitational mass with the most rapidly rotating model of neutron star in Table 3 of 
[jj; namely M = 1.864M . The radius used to describe the radial distances in the grid of our 
numerical integrations of z is r e = 10.755 km, which corresponds to the equatorial radius of the 
same rotating neutron star model. At the same diagram we have also plotted the relative difference 
between the exact metric component gfP and the one computed from transforming (r, 9) to (z, p) 
coordinates (lower thin solid curve). It is clear that the errors are at most of the order of ~ 10~ 6 in 
gu- Therefore any difference between a numerical metric and an analytical one should be of order 
at least 10 -6 to be considered true. The dashed vertical line marks the position on the z-axis that 
one reaches by integrating along the meridian of constant radius r = r$. The intersection of this 
line with the two curves depicts the systematic error induced in z (and correspondingly in gu) by 
replacing the dB/dx term in integration along the meridian by its numerical value. The errors 
at other locations along the z-axis are caused by numerical r-integration following the trapezoid 
rule. The deep wells in the plots are due to opposite signs of the errors accumulated along /x— and 
along r— integration in the region r > r$ that end up nullifying the difference between z^ and 
z^ at some point. Finally the vertical solid line on the leftmost part of the plot marks the polar 
surface location of the rotating neutron star model cited above in order to indicate the minimum 
z— value at which there is any physical meaning to look for numerical errors due to coordinate 
transformation. 
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a generic numerical metric, as long as the number of grid points in the vacuum region is 
sufficiently high (of the order of 300 x 300 for the the \i and the x coordinates in the vacuum 
region). To show this we shall appeal to the fastest rotating neutron star model with EOS 
FPS used in jl| (the model that shows up in the last line of the sequence with Mb = 2.105Mq 
of Table 3 of the corresponding paper) to estimate the magnitude of the error in computing 
the ^-coordinate. Of course in the case of a metric that is given in tabulated form we 
do not have a true value of the z coordinate to compare with, as in the simple example 
analyzed previously. On the other hand by a simple plot of B(xq,(jl) and (dB(x , fi)/dx)jq 
as a function of fi (where xq denotes the value of x-coordinate used for the first part of 
the integration path; namely the path along the meridian just outside the equator of the 
star, and the subscript n refers to the fact that the derivative is computed numerically), 
we conclude that both these functions are quite constant along the same meridian (B and 
(dB/dx)^ do not change by more than 0.1% and 25%, respectively, over the whole range 
of n), and thus the error in computing the first two integrals in ( l2TT) . even by the simple 
trapezoid rule is of the order of 

1 h 3 h 2 

(Az) M _ integration = -— /"(£) < — x max /"(//), (25) 
a 12 1U m 

where h denotes the step size of /i used in numerical integration (1/h is the total number of 
steps), while £ is some value of \i in the interval [0, 1]. We could form an upper value of this 
error by simplifying 1/12 to 1/10 and using the maximum value of f"(fi) instead of /"(£). 
It is easy to verify that this error in numerical integration over /i is of the order of only 
Az/M ~ 3 x 10~ 7 . This is systematically lower, or at most of the same order of magnitude, 
than the errors related to x-integration and numerical computation of dB/dx, that were 
estimated previously by means of the Schwarzschild example (c.f. Fig. [I]). Therefore, the 
plot of Fig. [1] summarizes quite well the overall numerical errors in computing the values 
of z and the corresponding errors induced in computing the g u component of a numerically 
constructed metric. 

We thus conclude that relative differences of metrics at the level of 10~ 6 and higher 
are true differences in metrics, and only these should be taken seriously into consideration 
when a proposed analytical metric is used as a faithful representation of a metric that is 
constructed numerically. The rest discrepancies between metrics could be easily attributed 
to inexact transformation of the coordinates. 
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V. CONCLUSIONS 



In this paper we have noted that when comparing numerical metrics with analytical met- 
rics, where the two metrics are expressed in different coordinates, one should be very careful 
in transforming coordinates. Especially since the Weyl-Papapetrou coordinates are usually 
used in vacuum Einstein stationary axisymmetric solutions, one should avoid the interior 
of a star in the integration path used to transform the initial quasi-isotropic coordinates, 
or else one results in erroneous comparison between metrics. Since the Cauchy-Riemann 
relations for the Weyl-Papapetrou coordinates do not hold in matter, integration of these 
relations along a path passing through the interior of the star results in a systematic shift 
of the z-coordinates, and consequently one ends up comparing two metrics at completely 
different points. 

By a thorough error analysis of the most stringent case (a maximally rotating neutron 
star model), we have concluded that when comparing a numerical metric with an analytical 
one, through faithful transformation of coordinates, any relative differences that exceed the 
level of 10 -6 should be considered real and not an artifact of the numerical transformation 
of coordinates, at least for metrics constructed on grids with grid size of at most that used 
by Stergioulas [ijj]. In Fig. [2] we have plotted once again the Figure 6 of 1], but now we have 
followed the path in the vacuum region just outside the star (the one described in Sec. IHIj) 
to integrate the coordinate z in order to transform the numerical metric components and 
compare them with the metric described by the solution of Manko et al. 8]. By direct 
inspection we find out that the matching between the two metrics is even better than what 
is inferred by Berti and Stergioulas PJ. The right transformation of the z-coordinate leads 
to about two orders of magnitude better matching between the two metrics than the one 
presented in even right at the surface of the star (leftmost part of the diagram). This 
renders the metric introduced by Manko et al. a very good candidate (even better than 
what it was first considered) to describe the space-time around a rotating neutron star. 

In a forthcoming paper we examine another similar candidate analytical metric (the one 
described in jzl, Q) to describe the space-time around any kind of neutron star, either 
rotating or not. Remember that the analytical metric used by Berti and Stergioulas had 
the disadvantage that it could not be adjusted to describe very slow rotating stars, since 
the corresponding metric could not be made to erase simultaneously both its quadrupole 
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FIG. 2: This is a replot of Figure 6 of [1] (upper curve) with somewhat altered dimensions along the 
horizontal axis (it's dimensionless here). The same difference between exactly the same two metrics 
is plotted for contrast (lower curve) when the transformation of coordinates is done according to 
our proposed method. The similarity between the analytic metric and the one at the exterior of the 
rapidly rotating neutron star model is manifestly much better than what was considered initially. 
The odd behavior of the lower curve near the surface of the star demonstrates a possible (probably 
accidental) crossing approach between the two metrics. The apparent oscillatory behavior of the 
rightmost part of the lower curve is due to the fact that the difference between the two metrics have 
approached the level of 10 -8 — 10~ 7 . This is actually the level of accuracy of the numerical method 
that have generated the numerical metric an thus spurious numerical information may have been 
introduced in the metric at this level. 

moment and its spin. 
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Appendix: The Schwarzschild metric in isotropic and Weyl-Papapetrou coordinates 

Here we write down the Schwarzschild metric in both isotropic and Weyl-Papapetrou 
coordinates in order to use it to measure the errors induced when (r, 9) — coordinates are 
numerically transformed to (p, z)— coordinates. Since in this example the exact metric is 
given as an analytical function in both sets of coordinates, we could infer what fraction 
of discrepancy between a numerical metric transformed into Weyl-Papapetrou coordinates 
and an analytical metric is caused by numerical errors in transforming the coordinates (see 
Sec. HE}. 

The spherically symmetric metric in Schwarzschild coordinates (t, f, 9, <fi) has the form: 
ds 2 = -(l- ™\ dt 2 +(l- 2 S] 1 dr 2 + r 2 (d9 2 + sin 2 9d<f) 2 ). (26) 



The same metric in isotropic coordinates (t, r, 9, <p) assumes the form: 

ds 2 = _ e 2^2 + e 2A ( rfr 2 + r 2 M 2 + r 2 ^2 q ^ ? (2?) 

where 

Jiv _ K 1 2r) „2\ _ / -1 , 1V1 \ 



The two radial coordinates, r, r are related to each other by 

From the form of the metric in isotropic coordinates (I28p one gets directly the B function: 



M 2 

B = l ~^T2- (30) 
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Finally, the same metric in Weyl-Papapetrou coordinates gets a complicated form, which 
could be obtained by writing down the inverse transformation of 

p = rsin^l-— J (31) 

z = rcos9(l + —) (32) 

and then constructing the metric components according to the recipe given in relations 
( |T4lll5|16p . Especially on the z-axis (9 = 0), Eq. ( 1321) could be easily solved with respect to 



2 

This function of z could then be used to obtain the exact value of all metric components as 
functions of ^-coordinate. For example the exact value of g u component, which we use in 
the paper to compare metrics, yields the following simple expression in Weyl-Papapetrou 

9„ M-z 

9tt = ~e 2u = — . 34 

y M + z v ' 
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